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(g) The disclosed technique and apparatus utTizes the different characteristic responses of a plurality of 
signals from a propagation or iriduction well logging device to determine the vertical and horizontal 
resistivity of formations surrounding a borehole. In a disclosed embodiment a method is set forth for 
determining the horizontal conductivity and the vertical conductivity of fonmations surrounding an 
earth borehole, including the following steps: a) deriving first and second fonmation conductivity 
values from measurements taken in the borehole; b) selecting a horizontal conductivity model value 
and a vertical conductivity model value ; c) computing error values from the differences between : (i) 
first and second composite conductivity values computed as a function of the horizontal and vertical 
conductivity model values, and (ii) the first and second derived formation conductivity values ; d) 
modifying the horizontal and vertical conductivity model values; e) repeating steps (c) and (d) until 
predetermined criteria of the error values are met ; and f) outputting modified horizontal and vertical 
conductivity model values. 
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HELD OF THE INVENTiON 

This invention relates to the field of earth borehole logging and. more particularly, to a method and appa- 
ratus for determining horizontal conductivity and vertical conductivity of earth formations. 

5 

BACKGROUND OF THE INVENTION 

It is recognized in the well logging art that formations surrounding an earth borehole can be anisotropic 
%vlth nagaid to conduction of electrical currents (see e.g. K.S. Kunz et ai., Some Effects of Formation Anisotropy 

10 on Resistivity Measurements In Boreholes, Geophysics, Vol. 23. No. 4, 1958). The phenomenon of electrical 
anisotropy is generally visualfa^d in one of fwo ways, or a combination thereof, as follows. 

In ntany sedimentary strata, electric current flows more easily in a direction parallel to the bedding planes 
than transversely to them. A reason for this anisotropy is that a great number of mineral crystals possess a fiat 
or elongated shape (e.g. mica and kaolin). At the time they were laid down, they naturally took an orientation 

IS parallel to the plane of sedimentation. The interstices in the formations are, therefore, generally parallel to the 
bedding plane, and the current is able to travel witfi facility along these interstices which often contain electri- 
caity conductive mineralized water. Such electrical anisotropy, sometimes called microscopic anisotropy. is ob- 
served mostly in shales. 

If a cylindrical sample is cut from a formation, parallel to the bedding planes, the resistivity of this sample 

20 measured with a current flowing along 'its axis is cafled the longitudinal (or horizontal) resistivity R|,. The inverse 
of Rf, is ti)e horizontal conductivity. Oh. If a simitar cylinder is cut perpendicular to fSne bedding planes, the re- 
sistivity measured witti a current flowing along its axis is called the transversal (or vertical) resistivi ty Ry. The 
inverse of Ry is the vertical conductivity, The anisotropy coefficient X, by definition, is equal to VR^(Or 

2s Voj/ct*,). Laboratory measurements have shown tiiat X may range from 1 to about 2.5 in different shales. 

Furtiiermore, the formations are often made up of a series of relatively thin beds having different lithdogic 
characteristics and, therefore, different resistivities (as. for example, sequences of thin shales and hanj 
streaks). In well logging systems the distances between the electrodes or antennas are great enough that the 
volume involved in a measurement may include several such tiiin beds. Since, in this situation, tiie current flows 

j0 more easily along the more conductive streaks tiian transversely to the series of beds, there is effective ani- 
sotropy. The effects on resistivity nteasurements of tiiis "macroscopic** anisotropy are cumulative with the ef- 
fects of the anisotropy due to the above-described mteroscopk; structure of the sediments. Reference can also 
be made to J.H. Moran et al., "Effects of Formation Anisotropy on Resistivity Logging Measurements. Geo- 
physics. Vol. 44, No. 7, 1 979, and to R. Chemali et al.. The Effect of Shale Anisotropy on FOcused Resistivity 

35 Devices", SPWLA Twenty-Eighth Annual Logging Symposium, 1 987. 

The determination of Rv as well as R^^ can be useful in various situations. For example, consider the case 
where the formation consists substantially of two types of material with resistivities and Rz with respective 
volume firactions « and 1-a. The effective horizontal and vertical resistivities Rf, and Ry&re given by 

40 ' (a/?2) + ((1-a)Rt) 

R, = (aR,) * {{^'a)R^ (2) 
If a is known, such as in a shale sequence where a gamma ray measurement or a spontaneous potential meas- 
urement has been used to provide the shale fraction. R^ and R2 can be determined from (1) and (2) if Rh and 
Rv are known. 

45 In situations where tiie borehole intersects the formations substantially perpendicular to the bedding 

planes, conventional induction and propagation well logging tools are sensitive almost exclusively to the hori- 
zontal components of the formation resistivity. When the borehole intersects the bedding planes at an angle. 
the tool readings contain an influence from the vertical resistivity as well as the horizontal. This is particularly 
true when the angle between borehole and the normal to the bedding planes becomes large, such as in direc- 

50 tional or horizontal drilling where angles near 90*^ are commonly encountered. In tiiese cases, tiie influence of 
vertical resistivity can cause discrepancies between measurements in these welts and measurements taken 
of the same formation in nearby vertical wells, thereby preventing useful comparison of these measurements. 
In addition, since reservoir evaluation is typically based upon data from vertical wells, use of data from wells 
drilled at high angles may produce erroneous estimates of formation produdbility if proper account is not taken 

55 of the anisotropy effect. 

A number of techniques have been proposed for measuring formation anisotropy andfor vertical conduc- 
tivity, such as by providing transmitter and/or receiver coils that are perpend icular to the borehole axis in addition 
to coils having conventional orientations. Reference can be made, for example, to U.S. Patent No.s 4.302.722, 
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4,302723, and 4.980.643. 

Equipment and techniques that determine horizontal and vertical conductivity (or anisotropy) by employing 
special equipment dedicated spedf icaily to such purpose result in increased equipment cost and incieased 
logging time and/or cost It is among the objects of the present invention to provide an apparatus and technique 
5 for determining horizontal and vertical conductivcty (or anisotropy determinable therefrom) using measure- 
ments that are often available from conventional types of equipment utilized for logging earth boreholes during 
drilling or by wireline. 

SUMMARY OF TOE INVENTION 

10 

A form of the present invention utilizes the different characteristic responses of a plurality of signals from 
a propagation or induction well logging device to determine the vertical and horizontal resistivity of formations 
surrounding a borehole. In a disclosed embodiment of the invention, a method is set forth for determining the 
horizontal conductivity and the vertical conductivity of formations surrounding an earth borehole, comprising 
15 the steps of: a) deriving first and second formation conductivity values from measurements taken in the bore- 
hole; b) selecting a horizontal conductivity model value and a vertical conductivity model value; c) computing 
error values from the differences between: (i) first and second composite conductivity values computed as a 
function of said horizontal and vertical conductivity model values, and (ii) said first and second derived forma- 
tion conductivity values; d) modifying said horizontal and vertical conductivity model values; e) repeating steps 

29 (c) and (d) until predetermined criteria of said error values are met; and f) outputting modified horizontal and 
vertical conductivity model values. {In the present application, any references to determination or use of resis- 
tivity are intended to generically means conductivity as well, and vice versa. These quantities are reciprocals, 
and mention of one or the ottier herein is for convenience of description, and not intended in a limiting sense.l 
The derived first and second formation conductivity values are obtained from measurements which are affected 

25 differently by the vertical and horizontal conductivities of the formations. In an embodiment of the invention, 
the horizontal and vertical conductivily model values are modified as a function of the error values. In this em- 
bodiment, the error values comprise a first error value that depends on the difference between the first com- 
posite conductivity value and ttie first derived formation conductivity value, and the second error value depends 
on the difference between the second composite conductivity value and the second derived formation conduc- 

30 tivf ty value. 

In a further embodiment of the invention, there is disclosed a method for determining the horizontal con- 
ductivity and the vertical conductivity of earth formations surrounding a borehole, including the following steps: 
a) deriving first and second formation conductivity values from measurements taken in the borehole, the derived 
first and second conductivity values being obtained from measurements which are affected differently by the 

35 vertical and horizontal conductivities of the formations; b) deriving an angle related to the dip of said formations; 
c) computing horizontal and vertical conductivity values from the derived first and second conductivity values 
and the derived angle; and d) outputting the horizontal and vertical conductivity values. In a form of this em- 
bodiment, the step of computing horizontal and vertical conductivity values includes storing a look-up table of 
values, and selecting values from the look-up table based on the derived first and second conductivity values 

40 and the derived angle. 

Further features and advantages of the invention will become more readily apparent from tfie following de- 
tailed description when taken in conjunction with the accompanying drawings. 

BRIEF DESCRIPTION OF THE DRAWINGS 

45 

Fig. 1 is a diagram, partially in block form, of an apparatus that can be utilized in practicing an embodiment 
of the invention. 

Fig. 2 is simplified diagram, partially in block form, of an apparatus for measuring formation resistivity that 
can be utilized in an embodiment of ttie invention. 
50 Fig.s 3 and a are, respectively, graphs of relationships between measured attenuation and formatk)n re- 
sistivity, and measured phase shift and formation resistivity. 

Fig. 5 shows a model space of the formations, with magnetic dipoles in the x and z directions. 

Fig.s 6Aand 6B, placed one below another, illustrate a routine for programming a processor in accordance 
with an embodiment of the invention. 
55 Fig. 7 is a simplified diagram, partially in block form, of an induction logging apparatus that can be utilized 

in an embodiment of the invention. 

Fig. 8 is a flow diagram of a routine for programming a processor to generate and store a look*up table 
useful in an embodiment of the invention. 
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Fig. 9 is a flaw diagram of a routine for programming a processor to use the ioolc-up table in accordance 
with an embodiment of the invention. 

DETAILED DESCRIPTION 

5 

Referring to Fig. 1, there is Illustrated an embodiment of the invention in conjunction with a logging-while- 
drilling apparatus. A platform and derrici< 10 are positioned over a borehole 11 that is formed in the earth by 
rotary drilling. Adrill string 11 is suspended within the borehole and includes a drSI bit 15 at Its lower end. The 
drill string 12 and the drill bit 15 attached ttiereto are rotated by a rotating table 16 (energized by means not 

10 shown) which engages a kelly 17 at the upper end of the drill string. The drill string is suspended from a hook 
18 attached to a travelling block (not shown). The kelly is connected to the hook through a rotary swivel 19 
which permits rotation of the drill string relative to the hook. Alternatively, the drill string 12 and drill bit 15 may 
be rotated from surface by a "top drive" type of drilling rig. 

Draiing f luki or mud 26 is contained in a pit 27 in ttie earth. A pump 29 pumps tiie drilling fluid into the drill 

15 string via a port in tiie swivel 19 to flow downward through tiie center of drill string 12. as indicated by flow 
arrow 30. The drilling fluid exits tiie drill string via ports in the drill bit 1 5 and tiien droulates upward in tiie region 
between the outside of the drilf string and tiie periphery of tiie borehole, as indicated by flow arrows 32. The 
drilling fluid tiiereby lubricates tiie bit and carries formatton cuttings to tiie surface of tiie earth. The drilling 
fluid is returned to the pit 27 for recirculation. 

20 Mounted in conjunction with the drill bit 15 is a bottom hole assembly 100 that may include an optional di- 
rectional drilling assembly 400 witti a mud motor having a bent housing or an offset sub (not separately rep- 
resented). Mounted above optional subassembly 400 is a formation resistivity (or conductivity) measuring ap- 
paratus 200 which is described fiirtiier hereinbelow. A measurement and communications subassembly 500 
is also provided and includes means, known in the art, for measuring and/or computing the direction and in- 

25 dination of the bottom hole assembly, the rotational orientation of ttie bottom hole assembly ("tool face"), and, 
typically, other measurements which need not be summarized herein. The communications portion of the sub- 
assembly 500 may typically be a mud pulse telemetry system that includes a downhole transmitter for sending 
coded information, including measurement signals, uphole and surface instrumentation for detecting and de- 
coding ttie information. One type of mud pulse telemetry system has an acoustic transmitter which employs a 

30 device known as a "mud siren" which includes a slotted stator and a slotted rotor that rotates and repeatedly 
interrupts the flow of drilling fluid to establish a desired acoustic wave signal in ttie drilling fluW. Driving elec- 
tronics in subassembly 500 may typically include a suitable modulator, such as a phase shiftkeying (PSK) mod- 
ulator, which conventionally produces driving signals for application to the mud transmitter. These driving sig- 
nals can be used to apply appropriate modulation to the mud siren. The generated acoustic mud wave travels 

35 upward in the fluid through the center of the drill string at the speed of sound in the fluid. The acoustic wave 
is received at the surface of the earth by transducers represented by reference numeral 31. The transducers, 
which are, for example, piezoelectric transducers, convert the received acoustic signals to electronic signals. 
The output of the transducers 31 is coupled to tiie uphole receiving subsystem 90 tiiat is operative to demod- 
ulate tiie transmitted signals, which can then be coupled to a processor 85 and a recorder 45. The processor 

40 85 can be of any suitable type, for example a digital microprocessor, with associated memory, timing, input/out- 
put, and display functk)ns. 

Fig. 2 is a simplified representation of the formation resistivity (or conductivity) measuring subassembly 
200 of the Fig. 1 apparatus. This subassembly, known in the art, may be in accordance witii apparatus described 
in U.S. Patent No. 4,899,112, which is incorporated herein by reference. In a preferred embodiment, tfie ap- 

45 paratus includes spaced apart transmitting antennas T, and Tj mounted on a section of drill collar 210, and 
spaced apart receiving antennas Ri and R2 mounted between the transmitting antennas. The antennas are 
coils supported in insulating media in recesses in tiie drill collar. The downhole electronics, typically contained 
in the drill collar 210 and represented in Fig. 2 by the block 220, indude, inter alia, means for alternately en- 
ergizing the transmitters with a signal having a frequency in the range about 0.1 MHz to 10 MHz, and typically 

50 about 2 MHz, Electromagnetic energy is transmitted into the formation and propagates therein. Energy shed 
back into the borehole is measured at the receiver pair, in a so-called differential receiver arrangennent, to ob- 
tain a phase shift measurement and an attenuation measurement. These measurements are coupled to the 
subassembly 500 for transmisskjn to the surface of the earth where means are provided, such as in processor 
85, for determining, as a function of the phase shift, ttie resistivity (Rp,) of formations at a relatively shallow 

55 depth of investigation around the receivers, and for determining, as a function of the attenuation, the resistivity 
(Rad) of formations at a relatively deep depth of investigation around the receivers. If desired, ttie resistivity 
values can be determined downhole and transmitted uphole for recording, or stored downhole for subsequent 
relrieval uphole. (The downhole subsystem furttier includes processor electronics, including a microprocessor 

5 
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with associated memory, timing, and interface circuitry, as described In the referenced U.S. Patent No. 
4,899,1 1 2.] Also, the use of two transmitters provides borehole compensation, but measurements could be ob- 
tained using a single transmitter. 

Fig 3 represents the relationship between the measured attenuation and the attenuation-dependent resis- 
6 tivity, Rad, as described in the above-referenced U.S. Patent No. 4.899,112. Fig. 4 represents ttie relationship 
between the measured phase shift and the phase-shift-dependent resistivity, Rp^ as described in said refer- 
enced Patent The illustrated exemplary graphs are for an embodiment with the receivers mounted on a 6.5 
inch diameter drill collar, receiver spacings of 25 and 31 inches, and a transmitter frequency of 2 MHz. As de- 
scribed in the referenced Patent, the plot of Fig. 3 can be utilized, for example, by entering the measured at- 
10 te nuation (minus a factor due to geometrical spread! ng loss) to read a corresponding value of Rad as determined 
by the curve. The plot of Fig. 4 can be utilized in similar fashion to obtain a value of Rp, lirom a phase shift 
measurement The referenced Patent indicates that a stored look-up table or stored polynomial curves can be 
used to automatically obtain values of Rad and/or Rp,, and that correction for permittivity can also be imple- 
mented. In the Rg. 1 arrangement, the processor 95 can be programmed in the manner described in the raf- 
fs erenced patent to obtain these values (and/or their inverses, the conductivities aad and <rps), or they can be com- 
puted downhole. 

Consider the model space in the diagram of Fig. 5 in which the z direction is considered vertical and the 
x-y plane considered horizontal. Initially, assume that bedding planes are perpendicular to the z direction and 
let Oy and cti, be the conductivities in the vertical (z) direction and the horizontal direction (x-z plane), respec- 
20 lively. As described in Moran et al., "Effects of Formation Anisotrapy on Resistivity Logging Measurements". 
Geophysics. Vol. 44. No. 7. 1979. the logging device coil fields can be considered as the superposition of mag- 
netic dipotes 511 and 512 (Fig. 5) having respective nr^gnetic moments. The magnetic dipole 511 is oriented 
in the vertical direction and the magnetic dipole 512 is oriented in a direction in the horizontal plane (arbitrarily, 
the X direction, in this case). The following notation is introduced. 
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Then the Hertz vector and scalar potential for a magnetic dipole located atx = y = z = Oand oriented in the x 
direction is given by: 
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For a similar dipole oriented in the z direction: 
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The electric and magnetic fields are given by: 



(6) 
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From these expressions for the electric and magnetic fields, the voltages at the receivers can b& derived 
as being equal to 



30 



35 



40 



where is the direction along the axis of the tool, is the radius of the receiver coil and R is the number of 
turns on the receiver. For an induction or a propagation logging tool M-nF^jTt where Rj is the radius of the trans- 
mitter coi, T is the number of turns on the transmitter coil and I is the current. 

From equations (6) and (7) one can write the magnetic field at the receiver coil for a transmitter at the origin 
and the tool tilted at an angle 0 from z in the x-z plane. For simplicity, we can assume that the too! remains in 
the x-z plane. Then the components of the dipole moment of the transmitter are given by: 

= MsinG (8) 
Ml = /Wcose (9) 

We will define H^ as the nnagnetic field in the i direction due to the component of the source in the j direction. 
This gives: 



45 



Xs 



(10) 



50 



55 



4n 



— k^xz —y- 



(11) 



7 



EP0 527089 A2 





e^^\ 3X2 


4Jt 


r» 1 



(12) 



Since the transmitter and receiver afe aligned in the same direction, the voltage In the receiver coll is given 

by: 

« = (hMiRRi(cos&(H^ + H^J + sln8(H^, + H^J) (14) 

Once the voltage at each of the receiver coils due to each of the transmitter coils is determined, the total meas- 
urement can be determined by adding the voltages in the case of an induction tool, or by taking the complex 
ratio of the voltages in the case of a propagation tool. For example, for the propagation logging device of Fig. 
2, the absolute value of the voltage at each receiver can be obtained as the square root of the sum of squares 

20 Of the real and imaginary parts of the complex voltage (equation (14)], and the ratio of the absolute values pro- 
vides the attenuation, from which the attenuation-determined resistivity R«, can be obtained (e,g. Fig, 3 above). 
The phase for each receiver is obtained from the arc tangent of the ratio of the inraginary and real parts of the 
complex voltage, and the phase shift Is the difference in phase at the two receivers. The phase-shift- 
determined resistivity Rp, can then be obtained (e.g. Fig.s 3 and 4 above). 

25 in accordance with an embodiment of the invention, formation conductivity values are obtained from meas- 

urements, as described in conjunction with Fig.s 1-4. From Uie two measurements and <jp^ an initial guess 
is made for the formation horizontal and vertical model conductivities and u^. The value(s) of the initial guess 
is not critical, and both can be tal^en to be equal to one of the measured values. Since the attenuation is gen- 
erally less affected by the anisotropy, it can serve as a good starting point Since both horizontal and vertical 

30 conductivities are equal for a starting point, the apparent values of a^^ and corresponding to these values 
of formation conductivity are simply equal to this value regardless of tool design or relative angle. If this agrees 
well with the observed values and o^^ then the formation may be judged to be isotropic and no further cal- 
culations are performed. 

If the initial model values do not fit the observed values, the next step involves computing the four partial 
35 derivatives dajdoh. da^jdcy . dcpjdch, and d^pjda^ These can be computed either analytically by talcing the 
derivatives of the equations above, or more easily by computing the finite differences between the values of 
<iad and Sps at slightiy different values of and Once the derivatives are computed, a step size and direction 
is chosen such that the step in and Ov will be along the direction to minimize the errors in the calculated 
values of and Cp,. The direction should be proportional to: 

Due to the extreme nonlinearity of the calculated values c«, and Op, upon the formation model values and 
it is important tiial the step size not be too large. It has been found that a good rule of thumb is that the 
maximum ratio d<j/cT should be less than 0.1 . In addition, the new sum of squared errors AoJ,+Aoi, is evaluated 
to be sure that It is not larger than the error from the previous step, where 

A<i^=cTp,-<Tps and Aa«<=CTa<rCTad- If ttie error is larger, tiien the step size is cut In half and another attempt Is made. 
If the new approximation is closer than the old one, then tiie process is repeated until a combination of at, and 
S„ is found which closely reproduces tfie observed values of and Wp,, or until a predetermined number of 
iterations have been performed without success. 

Referring to Fig. 6. there is shown a flow diagram of a routine for programming a processor (e.g. 85, Fig. 
1) to implement an embodiment of the invention. (It will be understood that the processing to be described need 
not be performed at the wellsite. and that signals obtained at the wellsite can be processed at a remote location.] 
The block 605 represents the reading in of the actual values and ap, determined from measurements at a 
particular depth level, and of a value for 9. such as from measured inclination obtained from measurement sub- 
assembly 500 of Fig. 1 . and/or from local knowledge. The block 608 represents the initalizing of the horizontal 
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conductivity and vertical conductivity model values. In the present embodiment the Initial values for both Ch 
and Gv are selected as being the measured conductivity value a^i for the present depth level, although it will 
be understood that other inrtiaf values can be utilized. The block 610 is then entered, this block representing 
inltalizing of "previous error values", A^jo and Apso, at an arbitrarily high Initial value, for purposes which will 

5 become understood. Also, an iteration counter is Initialized at zero. Block 612 is next entered, and represents 
the incrementing of the iteration count The block 615 is then entered, this block representing the compulation 
of the attenuation-determined conductivity and the phase-shift-determined conductivity, as a function of the 
model conductivity values and the dip, using equatk>n (14) in the manner prevtously described. The bk)ck 617 
te then entered, this block representing the computation of the emr component values Aed and Aps as the dif- 

10 ferences between the derived conductivity values and the computed conductivity values, inquiry is then made 
(diamond 620) as to whether both error values are less than 1 percent of the respective derived conductivity 
values. If not inquiry is made (diamond 622) as to whether a predetermined maximum iteration count (for ex- 
ample, 200 iterations) has been exceeded. If so, an indication of "no solution" Is provided for this depth level 
(block 625), and the diamond 627 is entered. If not a determination is made (diamond 630) as to whether the 

15 sum of squares of the error values Is greater than the sum of squares of the previous fold") error values. [For 
the first pass this cannot be the case sinoe the old error values were previously initialized at a high value.] If 
not, the block 633 is entered, this block representing the computation of the partial derivatives used in equations 
(15). and. subsequently, the computation of the step direction for incrementing the model conductivity values 
in accordance with equations (15) (block 642). Inquiry is then made (diamond 648) as to whether the deter- 

20 mined step sizes are greater than 10 percent of the present conductivity model values. If the condition is met 
for eitiier step size, the step sizes are proportionally reduced (blocks 652 and 654). so that the test of diamond 
648 would now be met So-called "old" nKxJel conductivity values are set equal to the present model conductivity 
values, as represented by Mock 668. The "old" conductivity model values are then incremented (block 672). 
and the "old" error values are set to the present error values (block 674). The block 612 is then re-entered for 

25 the next iteration. 

Returning, now, to the "yes" branch of diamond 630, when the sum of squares of the error values is deemed 
to be larger than l>efore, the step size is cut in half (block 631 ) and the conductivity model values are set equal 
to the "old" model conductivity values plus the step. The block 612 is ttien re-entered. If the test of diamond 
630 is still not met the step size is further reduced until the conditbn is satisfied. 

30 The described routine continues until either the errors are suf f icientiy minimized (condition of diamond 620) 
or the maximum number of iterations is exceeded (condition of diamond 622). If the error size condition is met 
the values of a, and ah are read out for the present depth level (block 640), e.g. to disc storage and recorder 
45 (Fig. 1), and the decision diamond 627 is entered. If the maximum iteration count is exceeded, an indication 
of no solution is read out (block 625) and decision diamond 627 is entered. If the last depth level to be processed 

35 has not been reached, the depth level index is incremented (block 629) and block 605 is re-entered to begin 
the process for the next depth level. 

Referring to Fig. 7, there is shown an embodiment of the invention wherein the derived conductivity meas- 
urement values are obtained from an inductk>n logging device 710 suspended from a wireline cable 780. The 
wireline cable is coupled with above-ground equipment that may include, for example, depth indicator 702, a 

40 receiver subsystem 790, and processor 95 and recorder 45 as in Fig. 1 . The subsurface induction logging equip- 
ment Is of a well known type, and only a simplified representation is set forth, showing a sonde 711 and trans- 
mitter and receiver coils T and R wound on insulating media on tiie sonde. The downhole electronics are rep- 
resented by the block 720. In induction logging, one or more transmitter coils (one of which is shown) are energ- 
ized by an AC signal. The resultant Oscillating magnetic field causes induction of currents in the formations 

45 which are approximately proportional to the formation conductivity. These currents, in turn, cause a voltage to 
be induced in one or more receiver coils (one of which is shown). Typically, the transmitter coil is energized by 
a signal generator, and the receiver coil is coupled to a drcuit that includes phase sensitive detectors which 
respectively receive reference phase signals that are in phase with the current in the transmitter coil and in 
phase quadrature witfi the current in the transmitter coil. Using these phase reference signals, tfie phase sen- 

50 sitive detectors respectively generate output signals which are proportional to the in-phase or "resistive com- 
ponent" of the induced receiver signal and the phase quadrature or "reactive" component of the induced re- 
ceiver signal. These output signals are used to obtain conductivity values respectively designated or and ax- 
It is known that tiiese conductivities are affected differentiy by the vertical and horizontal conductivities of the 
formations. [Reference can also be made to J.H. Moran et al,. "Basic Theory Of Induction Logging And Appli- 

55 cation To Study Of Two Coil Sondes", Geophysics . Vol. 27, No. 6, 1962. and to the publications and patents 
cited above.) The conduclivities gr and ax are illustrated as tjeing outputs of circuitry 720 that are transmitted 
to the earth's surface via cable 780. Processing can be performed downhole or uphole. Also, it will be under- 
stood that various known techniques, for example techniques relating to improved mutual signal cancellation 
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and focusing, can be employed. 

Processing for the case of and ax conductivity vaiues obtained from an Induction logging device is similar 
to that set forth In the flow diagram of Fig, 6 for the embodiment utilizing conductivity values obtained from a 
propagation logging device. In block 605, the values read in will be gr, the conductivity obtained from an in- 

5 phase induction logging measurement, and ax, the conductivity obtained from a phase quadrature induction 
logging measurement Regarding the computation of and ox (block 615 and subsequent), the voltage at the 
receiver(s) is determined from equation (14). and ap Is then obtained from the real part, with ox being obtained 
from the imaginary part When there are more than one transmitter coll and more than one receiver col, the 
receiver voltages attributable to each transmlttar-receh^er pair are added. The real part of the sum of voltages 

10 provides or, and tiie Imaginary part of the sum of voltages provides 0x- It will be understood that appropriate 
computational adjustments, to match the manner In which the measurements are taken and/or processed, can 
readily be made. 

It will be recognized that alternatives exist to the iterative approach described In conjunction witii the flow 
diagram of Fig. 6. For example, with derived values for the dip angle and for two conductivities from measure- 
is ments that are affected drfferentty by the vertical and horizontal conductivities of the formations (e.g. a^a and 
Op,, or Or and ox)' one can utilize the relationships set forth to obtain Oh and Oy in different manner. This may 
be done, for example, by generating and using look-up tables that are based on the relationships set forth or 
by using graphs representing the look-up tables. Polynomial functions, or other functional forms, whose coef- 
ficients are chosen in accordance with the relationships set forth, can also be utilized. 

20 Referring to Fig. 8, there is shown a flow diagram of a routine with which a processor can be programmed 

(typically, before operation) to obtain and store ttie look-up table for obtaining horizontal and vertical conduc- 
tivity values from input information tiiat Includes dip angle and Iwo measured conductivity values (e.g. and 
Op,, in tiiis example). A value of dtp angle is initialized atone end of tiie range the table is to cover, as represented 
by the block 811. Values of vertical conductivity and horizontal conductivity are also initialized at one end of 

25 ttie range that the table is to cover, as represented by the blocks 815 and 816. respectively. Values of and 
Ops are then computed, using equation (14) in the manner previously described, as represented by the block 
825. The present values of 0. ou and Ov are then stored at a memory location determined by the computed 
values and Ops. and by 0, as represented by the block 830. Inquiry is tiien made (diamond 835) as to whether 
the last Oh or the table has t>een reached. If not, o^ is incremented (block 638), Uie block 825 is re-entered, 

30 and the kx)p 840 continues until all Oh have been considered. When this occurs, inquiry is made (diamond 845) 
as to whether the last Oy in the range of interest has been considered. If not, Oy is incremented (tMock 854) Oh 
is re-initialized (block 818), and the nested loops 860 and 840 continue until all pairs of values of Oh and ot, are 
considered for the present 0. When tiie inquiry of diamond 845 is answered In the affirmative, table values will 
have been stored with regard to all pairs (oh, Oy) for the present value of dip angle, 0. Inquiry is then made (di- 

35 amond 862) as to whether the last 9 to be considered has been reached. If not, block 81 5 Is re-entered, Oy and 
Oh are again initialized, and the nested loops 840, 860. and 880 are repeated as table values are stored for 
each conductivity pair at each 6. 

Referring to Fig. 9. there is shown a flow diagram of a routine for controlling a processor, for example proc- 
essor 95, loaded with the previously computed look-up table information, to determine horizontal conductivity 

40 and vertical conductivity for particular derived values of aad. <7ps« and 9. The block 915 represents the reading 
in of the derived values of Oad, Ops, and 0 for a particular depth level. The closest address location in the data 
base is tiien determined (block 920). An interpolation is then performed (block 925) from said closest point 
(address), with respect to adjacent points in the three dimensions of the address space (block 925) to locate 
values of oh and Oy. The vaiues are read out (block 930), and inquiry is made (diamond 940) as to whether the 

45 last depth level to be processed has been reached. If not, the next depth level is considered (block 950), and 
the t>lock 915 is re-entered for processing the next depth level. The procedure continues until the inquiry of 
diamond 940 is answered in the affirmative. 

The invention has been described with reference to particular preferred embodiments, but variations witiiin 
the spirit and scope of the invention will occur to those skilled in the art. For example, it will be understood that 

50 the derived measurements, which respectively have different dependence on the formation horizontal and vert- 
ical conductivities, can be from other types of logging devices, such as devices employing electrodes. Also, it 
will be understood that the measurements can be from different devices, for example one measurement from 
a propagation logging device and another measurement from a logging device that utilizes electrodes. Further, 
tiie respective measurements can be taken at different frequencies, different spacings. or with otiier parameter 

55 variations. It should also be noted that while two conductivity measurements are utilized in the illustrated em- 
bodiment, more than two such measurements can be utilized at a particular depth level, if desired. Also, the 
principles of the present invention are applicable where anisotropy exists in three dimensions, for example with 
different conductivities in the x.y directions in the horizontal plane. Finally, it will be understood that reference 
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to honzontat and vertical components is for convenience, and is intended to generically include processing to 
obtain orthogonal conductivity components. 



5 Claims 

1. A metiiod for determining the horizontal conductivity and the vertical conductivity of earth formations sur- 
rounding a borehole, the method including the step 

(a) of deriving first and second formation conductivity values from measurements taken in the borehole, 
10 and charactarized by the steps of: 

b) selecting a horizontal conductivity model value and a vertical conductivity model value; 

c) computing error values from the differences between: (i) first and second composite conductivity val- 
ues computed as a function of said horizontal and vertical conductivity model values, and (i!) said first 
and second derived formation conductivity values; 

15 d) modifying said horizontal and vertical conductivity model values; 

e) repeating steps (c) and (d) until predetermined criteria of said error values are met; and 

f) outputting modified horizontal and vertical conductivity model values. 

2. The method as defined by claim 1, wherein said derived first and second formation conductivity values 
20 are firom respectively different measurements taken at a given depth level in the t>orehole. 

3. The method as defined by claim 2, wherein said derived first and second formation conductivity values 
are obtained from measurements which are affected differently by the vertical and horizontal conductiv- 
ities of the formations. 

25 

4. The method as defined by daim 2 wherein said derived first and second formation conductivity values 
are respectively obtained from attenuation and phase measurements. 



JO 



50 



5. The method as defined by claim 2, wherein said first conductivity value is derived from an in-phase in- 
ductton logging measurement in ttie borehole and said second conductivity value is derived from a phase 
quadrature inductton logging measurement in the tsorehole. 



6. The method as defined by claim 5, further comprising deriving a dip angle of the boretiole with respect to 
the formation; and wherein said step (c) includes computing an in-phase conductivity value and a phase 
quadrature conductivity value as functions of said horizontal and vertical conductivity model values and 
said dip angle; and computing an error as a function of the derived first and second conductivity values 
and the computed conductivtty values. 

7. The method as defined by claims 1 , 2. 3 or 4, wherein said horizontal and vertical conductivity model val- 
ues are modified as a function of said error values. 

40 

8. The method as defined by claims 1, 2, 3, 4 or 7, further comprising deriving a formation dip angle, and 
wherein said first and second composite conductivities are computed as a function of said horizontal and 
vertical conductivity model values and sakl dip angle. 

^5 9. The method as defined by any of claims 1 -4, wherein at least said steps (c) through (e) are performed by 
a machine processor. 

10. The method as defined by daim 2, wherein said horizontal and vertical conductivity models are modified 
as a function of said error values, and wherein said error values comprise a first error value that depends 
on the difference between said first composite conductivity value and said first derived formation conduc- 
tivity value, and said second error value depends on the difference between said second composite con- 
ductivity value and said second derived formation conductivity value. 



1 1 . The method as defined by daim 1 0. wherein said predetermined criteria of said error values requires both 
said first and second error values to be less than prescribed maxima. 

55 

12. The method as defined by any of daims 1-4. wherein said steps (c) through (e) are repeated until said 
errorvalues are less than prescribed maxima or until said steps have been repeated a predetermined num- 
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ber of times. 

13. The method as defined byanyofclaims 1-4, further comprising repeating steps (a) through (f) for further 
depth levels in the borehole, and producing a recording of the outputted conductivity model values. 

14. A method for determining the horizontal conductivity and the vertical conductivity of earth formations sur- 
rounding a t)orehole. the method including the step 

(a) of deriving first and second formation conductivity values from measuren^nts taken in the borehole, 
said derived first and second conductivity values being obtained from measurements which are affected 
differently by the vertical and horizontal conductivities of the formations; and characterized by the steps 
of: 

b) deriving an angle related to the dip of said formations; 

c) computing horizontal and vertical conductivity values from said derived first and second conductivity 
values and said derived angle; and 

d) outputting said horizontal and vertical conductivity values. 

15. The method as defined by daim 14, wherein said step of computing horizontal and vertical conductivity 
values includes storing a look-up table of values, and selecting values from said look-up table based on 
said derived first and second conductivity values and said derived angle. 

16. The method as defined by daim 14 or 1 5, wherein said derived first and second conductivity values are 
respectively an attenuation-dependent conductivity vahje and a phase-shift-dependent conductivity value 
obtained from a propagation logging device. 

17. The method as defined by daim 14 or 15, wherein said derived first and second conductivity values are 
respectively an in-phase conductivity value and a phase-quadrature conductivity value obtained from an 
induction logging device. 

18. The method as defined by daim 14 or 15, wherein said steps (a) through (d) are repeated for a number 
of depth levels in the borehole to produce a recording of outputted conductivity values. 

19. Apparatus for determining the horizontal conductivity and the vertical conductivity of earth formatk)ns sur- 
rounding a borehole, including means for deriving first and second formation conductivity values from 
measurements taken in the t>orehole, and characterized by: 

means for selecting a horizontal conductivity model value and a vertical conductivity model value; 

means for computing error values from the differences between: (i) first and second composite con- 
ductivity values computed as a function of said horizontal and vertical conductivity model values, and (ii) 
said first and second derived formation conductivity values; 

means for modifying said horizontal and vertical conductivity model values; 

means for repeating implementation by said computing means and said modifying means until pre- 
determined criteria of said error values are met; and 

means for outputting modified horizontal and vertical conductivity model values. 

20. Apparatus as defined by daim 19, wherein said derived first and second formation conductivity values 
are from respectively different measurements taken at a given depth level in the borehole. 

21. Apparatus as defined by daim 20, wherein said derived first and second formation conductivity values 
are obtained from measurements which are affected differently by the vertical and horizontal conductiv- 
ities of the formations. 

22. Apparatus as defined by claims 1 9. 20 or 21 , further comprising deriving a formation dip angle, and where- 
in said first and second composite conductivities are computed as a function of said horizontal and vertical 
conductivity model values and said dip angle. 
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Fig.l 




13 



EP 0 527 089 A2 



Fig. 2 
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Fig. 4 
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Fig. 6A 
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Fig, 8 
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